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Abstract 

The effective potential for the Standard Model Higgs field allows two quasi-degenerate vacua; 
one is our vacuum at the electroweak scale, while the other is at a much higher scale. The latter 
minimum may be at a scale much smaller than the Planck scale, if the potential is lifted by new 
physics. This gives rise to a possibility of domain wall formation after inflation. If the high- 
scale minimum is a local minimum, domain walls are unstable and disappear through violent 
annihilation processes, producing a significant amount of gravitational waves. We estimate the 
amount of gravitational waves produced from unstable domain walls in the Higgs potential and 
discuss detectability with future experiments. 
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I. INTRODUCTION 


The Standard Model (SM) of particle physics has been extremely successful, and the 
Higgs boson with a mass of 125 GeV discovered at the LHC [1, 2] completed the last 
missing piece of the SM. So far there is no experimental hint for physics beyond the SM, 
and it may be that the Standard Model is valid up to very high energy scales beyond the 
reach of the current collider experiments. 

In the SM framework, the measured values of the Higgs boson mass and top quark 
mass imply that the electroweak (EW) vacuum is likely nretastable. This is because the 
Higgs self coupling becomes negative at some high energy scale as a result of the fact that 
its RGE (renormalization group equation) evolution is dominated by the top Yukawa 
coupling [3, 4], While such effective potential is acceptable as long as our vacuum is 
sufficiently long-lived, it may signal that new physics appears around that scale and lift 
the potential. 1 For example, higher dimensional operators may lift the effective potential 
so as to make these two vacua degenerate in energy, or even make the EW vacuum stable. 

The existence of two quasi-degenerate minima in the Higgs potential has interesting 
cosmological implications. During inflation, either of the two vacua is randomly selected 
in each patch of the Universe, if the Higgs field acquires quantum fluctuations large enough 
to overcome the potential barrier between the two minima. Then, domain walls are formed 
after inflation. Domain walls are a sheet-like topological defect [5] , and they are stable if 
the two vacua are exactly degenerate. Stable domain walls, however, are a cosmological 
catastrophe as they eventually dominate the Universe, generating unacceptably large 
inhomogeneities. 2 If there is an energy difference (bias) between the two vacua, domain 
walls are unstable and they eventually annihilate [11-14]. The EW vacua is realized in 

1 The scale sensitively depends on the top quark mass, and it will be around or beyond the Planck scale 
for the top quark mass about 171 GeV, the lower side of the experimental range [3, 4]. 

2 If the high-scale minimum is around or beyond the Planck scale, an eternal topological Higgs inflation 
may take place avoiding the cosmological disaster [6]. See Refs. [7, 8] for the original works of topological 
inflation. The topological Higgs inflation may provide a dynamical explanation for the multiple-point 
criticality principle [9]. See also Ref. [10] for the related topics. 
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the whole Universe if it is energetically preferred. Interestingly, a significant amount of 
gravitational waves is emitted in the violent annihilation processes [15-18]. As we shall 
see shortly, such gravitational waves are within the sensitivity reach of future experiments 
such as advanced-LIGO [19], KAGRA [20, 21], ET [22], LISA [23] and DECIGO [24], if 
the domain walls are sufficiently long-lived. Thus, gravitational waves can be a probe of 
another vacuum far beyond the EW scale. 3 

In this letter, we study the gravitational waves generated by collapsing domain walls 
in the Higgs potential with the quasi-degenerate vacua. In Sec. II, we introduce the Higgs 
potential having the false vacuum at high energy scales and discuss the possibility of the 
domain wall formation. In Sec. Ill, we calculate the gravitational wave abundance and 
discuss its detectability with future experiments. Sec. IV is devoted to discussion and 
conclusions. 

II. FALSE VACUUM IN HIGGS POTENTIAL 

Let us consider the following Higgs potential lifted by new physics at some high energy 
scale, 

V„ = jAfoV + |[[. (1) 

where ip is the Standard Model Higgs scalar field, A(^) is a scale-dependent self coupling 
constant and A is a cutoff scale for the dimension six operator. We have neglected the 
quadratic term of order the electroweak scale as we are interested in the behavior of the 
Higgs potential at high energy. Here we simply substitute the Higgs field value for the 
renormalization scale, A (p) = A {ip). 

The renormalization-group-improved effective potential including one-loop and two- 
loop corrections in Landau gauge is given in [3] , and the gauge-dependence of the effective 
potential was examined in [4], More recently, a treatment for higher dimensional operators 
in the Higgs potential was studied in [27]. For our order of magnitude estimate, however, 

3 The gravitational waves from domain wall annihilation can also be a probe of the SUSY breaking 
scale [25] or a thermal inflation scenario [26]. 
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the following crude approximation is sufficient. We leave further refinement of our analysis 
for future work, but we believe our main results will not be qualitatively changed. 

The scale dependence of the Higgs self coupling A (p) is governed by the RGE, 

( 2 ) 

where t = In (n/M t ) and 

f3\ = 24A 2 + 12A yf — Qyf — 3A (g r2 + 3 g 2 ) + -[2 g 4 + ( g ' 2 + g 2 y] (3) 

up to one loop order. Here y t is the top Yukawa coupling, M t is the top quark mass, 
g' and g are respectively the U(1)y and the SU(2) L gauge coupling constants. The top 
Yukawa coupling gives a dominant contribution to the RGE up to a very high energy 
scale where the U(l)y gauge coupling becomes large. As a result, the Higgs self coupling 
turns to negative, and the effective potential becomes negative at an intermediate scale 
in the absence of higher dimensional operators. 

The effective potential can be lifted by higher dimensional operators such as the last 
term in the right-hand-side in Eq. (I). 4 In this case there are two potential minima; one 
is at the EW scale <p = t> EW , an d the other at a much higher scale ip = ipf. Depending 
on the size of the higher dimensional operator, the high-scale minimum can be a local or 
global minimum. In particular, our main interest lies in the case when the two minima 
are quasi-degenerate and the EW vacuum is slightly energetically preferred: 

V H {v ew) ~ V H (<Pf) ~ o, (4) 

Vr{v ew) < Vff (<£>/)■ (5) 

If the Higgs field acquires a sufficiently large quantum fluctuations during inflation, both 
vacua may be populated in different patches of the Universe, leading to domain wall 
formation after inflation. In a later Universe, the EW vacuum will be selected after 
domain wall annihilation. 

4 Alternatively, one may lift the potential by introducing e.g. an additional singlet scalar field which 
gives a positive contribution to the beta function (3) [29-31]. Our results hold in this case too, without 
significant modifications. 
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FIG. 1: The effective potential for the Higgs field, Vh{< p). We have taken M t = 173.28 (solid 
red), 173.29 (dashed green) and 173.30 GeV (dotted blue), mn = 125.6 GeV, and A = 10 10 
GeV. Note that the values of A and ipf are for illustration purposes only. See the text for details. 

In Fig. 1 we show the Higgs potential calculated in the present set-up for illustration 
purposes. We have taken the Higgs boson mass rn h = 125.6 GeV, M t = 173.28 (solid 
(red)), 173.29 (dashed (green)), 173.30 GeV(dotted (blue)) and A = 10 10 GeV. One can see 
that the position of the high-scale minimum, tp — iff, is comparable to, but parametrically 
smaller than A because of the loop suppression factor in the RGE. The cut-off scale as 
well as (fif increase for smaller values of M t . For notational simplicity, let Vf and V max 
denote the potential energy at the false vacuum and the local maximum, respectively, as 
shown in the figure for the middle line. 

Note that the precise value of A obtained in Ref. [4] is about two orders of magnitude 
larger and it is about 10 12 GeV for M t = 173.3 GeV. Accordingly, iff will be larger by a 
similar amount. We emphasize here that our analysis in the next section does not depend 
on the detailed RGE evolution, because we express all the relevant quantities in terms of 
Pf, Vmax and Vf. One should simply use the result of e.g. Ref. [4] when one relates the 
value of iff to the top quark mass and the Higgs boson mass. 
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III. GRAVITATIONAL WAVES FROM COLLAPSING DOMAIN WALLS 


As we have seen in the previous section, the Higgs potential allows two quasi-degenerate 
minima, especially if the potential is lifted by higher dimensional operators. This gives rise 
to a possibility of domain wall formation after inflation. A domain wall is characterized 
by its tension, a, which is roughly estimated to be 

^ ~ + ^ max ) w ~ ^ l<Pf, ( 6 ) 

where 14 aax is the height of the potential barrier between two minima, w ~ (pf/V m kx is the 
width of the domain wall, and we fix it to minimize the tension in the second equality. 
In order to avoid the cosmological domain wall problem, the energy bias between the two 
vacua is necessary to make domain walls unstable. In the presence of the bias, domain 
walls start to collapse when the energy density of the domain walls become comparable 
to the bias energy density. As is confirmed by numerical simulations, the evolution of the 
domain wall network exhibits a scaling behavior [32-35] and the energy density of the 
domain wall is roughly given by 

Pdw ~ cr H. (7) 


Then, the Hubble parameter at the domain wall decay is 


H V f 


Vf 


cr 


Vina 


( 8 ) 


In order not to generate unacceptably large inhomogeneities, domain walls must decay 
before they dominate the Universe, which places a lower bound on the energy bias; 

ffd “ > Hd ° m ^ vd > {mp) • (9) 

where M P ~ 2.4 x 10 18 GeV is the reduced Planck mass, and H dom ~ cr/Mj, is the Hubble 
parameter when domain walls would start to dominate the energy density of the Universe 
if there were not for the bias. 

The domain wall collapse is a violent processes, and some part of the energy stored in 
the domain walls is converted to gravitational waves. The spectrum of the gravitational 
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waves is expected to be peaked at a frequency corresponding to the Hubble scale at the 
decay, as it is the typical curvature scale of the domain wall system. This was confirmed 


by detailed numerical calculations [18], and the density parameter of the gravitational 
waves at the peak frequency at the time of domain wall collapse is given by 



( 10 ) 


where G is the Newton’s gravitational constant and and A are numerical factors 
characterizing respectively the efficiency of the gravitational wave emission and the area 
of the domain walls. They are determined by the numerical calculations to be e gw ~ 0.7 
and A — 0.8 [18]. Then, the present time density parameter of the gravitational waves at 
the peak frequency is obtained as 



( 11 ) 


where g * is the relativistic degrees of freedom at the domain wall decay and 7 is the 
dilution factor after the domain wall decay, and it is given by 



for H dec < H r 


1 


( 12 ) 


(H R /H dec y/ 3 for H dec -> Hr 


Here Hr is the Hubble parameter at the reheating. The peak frequency corresponds to 
the Hubble parameter at domain wall decay, which is red-shifted by the cosmic expansion 
until today, 



(13) 


where 



-^dec for //dec < Hr 
Tr for // dec > Hr 


(14) 


and Xdec and Tr are the cosmic temperature at H = H dec and Hr, respectively. 

Now let us turn to the domain walls in the SM Higgs potential lifted by new physics. 
We focus on the case in which the high-scale minimum is a false vacuum and it is quasi¬ 
degenerate with the EW vacuum. The position of the false vacuum is at intermediate 
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energy scales, 10 s GeV < <p/ < 10 12 GeV, depending on the values of the top quark mass, 
the strong gauge coupling, and the Higgs boson mass [3, 4], The height of the potential 
barrier is determined by solving the RGE, which is roughly V max rs-/ 10 A ip A . The bias 
energy density is treated as a free parameter which is adjusted by tuning A. 

In the case of the Higgs domain walls, there are generically finite temperature correc¬ 
tions to the Higgs potential, which give an extra contribution to the (time-dependent) 
energy bias. In the following we derive a condition for thermal effects to have a negligible 
impact on the domain wall dynamics. The condition can be relaxed in certain situations, 
and we shall give concrete examples later. 

Throughout this letter, we focus on the case in which the position of the false vacuum 
is always much larger than the cosmic temperature, ipf T. Then the thermal mass 
correction to the effective potential is negligibly small at p — ipf. Even for p /§> T, 
however, there is a logarithmic correction arising from the free energy of thermal plasma, 
the so called thermal log potential [36] , which is roughly given by 

V T {ip) = aT i (15) 


where a is a numerical constant of 0(0.1). In addition, there are background thermal 
plasma in the EW vacuum, while many of the SM particles are non-relativistic in the 
false vacuum because of ipj T. These thermal effects are considered to generate an 
extra energy bias of order T 4 . 

After reheating, the thermal energy bias ~ T 4 decreases faster than the energy density 
of domain walls in the scaling regime. Before reheating, both evolve in the same way since 
the dilute plasma energy density evolves as T 4 ~ T^HMp oc pd w (see Eq. (7)) in the case 
of usual perturbative decay of the inflaton. Thermal effects do not induce the domain 
wall annihilation if pd w ~ > T 4 . For a given reheating temperature, this gives a lower 

bound on the field value at the false vacuum: 


Tr < 3 x 10 8 GeV 


<Pf \ 

10 12 GeV ) 


3/2 



(16) 


where the second inequality represents the condition for avoiding the domain wall domi¬ 
nation (9). The condition (16) implies that the domain walls must annihilate before the 



reheating. 

For the reheating temperature satisfying the condition (16), the thermal plasma energy 
is always much smaller than F max , and so, the use of the tension given in Eq. (6) is 
justified. 5 If the condition (16) is violated, domain walls will disappear soon after inflation 
and the resultant gravitational wave signal will be too weak to be detected by future 
experiments. Also the peak frequency tends to be extremely high. 

We would like to emphasize that the bound (16) should be taken with care, because 
the thermal history after inflation is unknown. Indeed, it is possible to consider some 
complicated (and contrived) thermal history where thermal effects on the domain wall 
dynamics are negligible even if (16) is not satisfied. For instance, the inflaton may domi¬ 
nantly decay into hidden sector particles, and the SM sector is reheated when the coupling 
between the SM and hidden sectors freezes-in well after the domain wall annihilation. 6 
We shall return to this issue in Sec. IV. 

In Fig. 2, we show the contours of the gravitational wave density parameter at the 
peak frequency, flcw^lpeak (dashed (red) lines), as well as the peak frequency, / pea k 
(dotted (blue) lines), in the plane of (iff, (Ify/fynax) 1 / 4 ). We have set Tr — 3 x 10 s GeV 
in Fig. 2(a) and Tr = 10 4 GeV in Fig. 2(b). The shaded (magenta) lower-right triangle 
region is excluded by the domain wall domination (cf. Eq. (9)) and the solid cyan line is 
the lower bound on iff to avoid the early domain wall decay due to the thermal effects 
(cf. Eq. (16)). As mentioned earlier, thermal effects are negligible in the region right to 
the solid cyan line, but it does not necessarily preclude the region left to the solid cyan 
line because of large uncertainties of thermal history. The thick green and yellow lines 
represent the sensitivity curves of advanced-LIGO and ET respectively in Fig. 2(a), and 
LISA and DECIGO in Fig. (2(b)), and the shaded region below the curves will be probed 

5 The evolution of the domain wall network might deviate from the scaling regime in the presence of 
background plasma. The existence of plasma could also change the effective area of domain walls A, 
which results in a shift of the peak frequency for the fixed tension and energy bias. We however expect 
that such thermal effect on the domain wall dynamics is small as long as thermal energy is much smaller 
than the energy stored in the domain walls, as we assume in the text. 

6 Our estimate (10) remains unchanged in this case. 
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(a) T r = 3 x 10 s GeV 



(b) T r = 10 4 GeV 

FIG. 2: The contours of ^GW^Ipeak (dashed red) and / pea k (dotted blue) in the plane of 
<Pf~(Vf/ Fmax) 1 ^ 4 - We have set Tr = 3 X 10 8 GeV and Tr = 10 4 GeV in the upper and lower 
panel, respectively. The dash-dotted black line represents the border above (below) which the 
reheating takes place after (before) the wall decay (16). The lower right shaded (magenta) 
region is excluded by the domain wall domination (9) and the solid cyan line corresponds to the 
lower bound on (ff due to the early domain wall decay by the thermal effect. The thick green 
and yellow lines represent the sensitivity curves for advanced-LIGO (LISA) and ET (DECIGO) 
respectively in the upper (lower) panel and the shaded region below the curves will be probed 
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by each experiment. 





FIG. 3: The typical spectrum of the gravitational waves is shown by the solid (red) lines. We 
have taken pf = 2x 10 9 GeV and (V//Vmax) 1//4 = 5x 10 -5 for the left line and iff = 2x 10 12 GeV 
and (Vf/Vmax) 1 / 4 = 10 -3 for the right line. 

by each experiments. The sensitivity curve of KAGRA is expected to be similar to that of 
advanced-LIGO. Note that it is difficult to generate gravitational waves within the reach 
of pulsar timing observations such as IPTA [37] and SKA [38] which are sensitive to much 
lower frequencies (~ 10 -9 Hz). It would require the tension of the domain walls to be 
close to (10 5 GeV) 3 , which cannot be realized within the SM framework. 

We show in Fig. 3 the gravitational wave spectrum for certain parameters (pj = 
2 x 10 9 GeV and (Vf/V max ) 1 ^ = 5 x 10 “ 5 for the left line and ipf = 2 x 10 12 GeV and 
{Vf/V max ) 1 / 4 = 10 ~ 3 for the right line). According to the precise numerical calculations 
[18], the gravitational wave spectrum scales as (/// pea k ) 3 for / < f peak and (///peak )^ 1 
for / > / peak . 

IV. DISCUSSION AND CONCLUSIONS 

In the previous section we have mentioned that the condition (16) can be relaxed in 
certain situations. Here we give concrete examples that relax the condition (16) so as to 


11 







make the wider parameter available. As we shall see below, this requires a more involved 
thermal history of the Universe. 

Suppose that the inflaton dominantly decays into a hidden sector, which is decoupled 
from the SM sector. This significantly suppresses thermal corrections to the Higgs po¬ 
tential, and the lower bound on tpf (16) is no longer applied.' The SM particles can be 
thermally populated through various processes much later than the domain wall annihila¬ 
tion. For instance, the hidden sector may contain a U(l) gauge group, which has a small 
kinetic mixing with the SM U(l)y. For a sufficiently small kinetic mixing, the SM par¬ 
ticles are thermalized through the kinetic mixing well after the domain wall annihilation 
but before nucleosynthesis, unless the hidden photon is extremely light. If the hidden 
photon is massless, we need to introduce hidden matter fields (e.g. a Dirac fermion) 
charged under the hidden U(l) gauge symmetry. The hidden matter fields behave as 
mini-charged particles, connecting the hidden sector and the SM sector. In this case, the 
massless hidden photon contributes to the effective neutrino species, N e $, whose precise 
value depends on the mass of the hidden matter fields [39]. Similarly, if one introduces a 
gauge singlet scalar field, the SM sector can be produced from the hidden sector through 
the Higgs portal coupling. In this case, the thermalization of the SM sector can be delayed 
for a sufficiently small Higgs portal coupling. It is worth noting that there is no entropy 
production when the SM sector is thermalized, and our analysis in the previous section 
can be applied without any modification even in the parameter space left to the solid cyan 
line in Fig. 2. 

Alternatively, the SM sector can be reheated, if one introduces another held (e.g. a 
modulus held) which dominates the Universe for a short period and decays into the SM 
particles after the inhaton decay. In this case one has to take account of its thermal 
effects, the shift of the peak frequency, as well as an extra dilution due to the modified 

' Some amount of the SM plasma is necessarily generated by the Higgs domain walls, because they 
continuously annihilate in the scaling regime. We have numerically confirmed that the energy density 
of the SM plasma induced by the Higgs domain walls in the scaling regime is always smaller than that 
of domain walls, and therefore, it does not induce the domain wall annihilation. 
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thermal history. 

In the SM framework, the Higgs self coupling turns to negative, and the effective po¬ 
tential becomes negative at a high energy scale, based on the perturbative RGE analysis. 
The scale sensitively depends on the values of the top quark mass, but it is at an inter¬ 
mediate scale for the top quark mass, M t ~ 173 GeV. In such a case, the EW vacuum is 
metastable, which is acceptable as long as it is sufficiently long-lived. On the other hand, 
the negative effective potential may signal that new physics appears around that scale 
and lift the potential, creating a local minimum at an intermediate scale. In a limiting 
case, the two vacua, one at the EW scale and the other at an intermediate scale, are 
quasi-degenerate in energy. This requires a certain amount of fine-tuning of the param¬ 
eters. Once the fine-tuning of the parameters is realized, it gives rise to a possibility of 
domain wall formation. 

For domain walls to be formed after inflation, both vacua must be populated with more 
or less equal probability during inflation. Domain walls are not formed if one of the two 
vacua is preferred over the other. This is the case if the Higgs potential is significantly 
modified through its non-minimal coupling to gravity, or if the quantum fluctuations of 
the Higgs field is too small to overcome the barrier separating the two vacua. We have 
assumed that the Higgs field is minimally coupled to gravity as well as the inflaton sector 
so that the potential is not significantly modified, and that the quantum fluctuations of 
the Higgs field is sufficiently large to overcome the barrier between the two vacua. 

The energy of domain walls is partially converted to gravitational waves through the 
violent annihilation processes. The gravitational waves are expected to be peaked at a 
frequency corresponding to the Hubble horizon size at the domain-wall annihilation. The 
gravitational-wave spectrum looks like that from a phase transition. In the latter case, 
the gravitational wave spectrum has a peak corresponding to the typical size of bubbles, 
a few orders of magnitude smaller than the Hubble horizon at the phase transition. 

In order to generate gravitational waves within the reach of future experiments, the 
bias energy density must be so small that domain walls annihilate when they are about 
to dominate the Universe. While there is no compelling reason for the quasi-degeneracy 
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of the two vacua, it may be due to the multiple-point principle [9] or an important 
cosmological role (e.g. dark matter, baryogenesis) of the decay products of domain walls. 
The violent domain wall annihilation processes produce not only gravitational waves but 
also a large amount of the SM Higgs bosons which soon decay into quarks, leptons and 
gauge bosons. If the SM Higgs boson is coupled to some heavy degrees of freedom such 
as right-handed neutrinos, a B-L Higgs held, and a singlet scalar, those heavy particles 
can be produced through non-perturbative processes, and they may contribute to dark 
matter, or baryogenesis, etc. [40]. 

In this letter we have investigated the domain wall formation in the Standard Model 
Higgs potential lifted by new physics and gravitational wave production due to the domain 
wall decay. We have shown that the gravitational waves can be within the reach of the 
future experiments such as advanced LIGO, KAGRA, ET, LISA and DECIGO, if the 
domain walls decay when their energy density is sizable. In this way, the direct detection 
experiments of the gravitational wave can be a powerful probe of another vacuum far 
beyond the EW scale. 
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